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ABSTRACT 

A simple model is constructed to evaluate the change of star formation rate of 
a disk galaxy due to environmental effects in clusters of galaxies. Three effects, 
(1) tidal force from the potential well of the cluster, (2) increase of external 
pressure when the galaxy plows into the intracluster medium, (3) high-speed 
encounters between galaxies, are investigated. General analysis indicates that 
the star formation rate increases significantly when the pressure of molecular 
clouds rises above ~ 3 x 10^ cm~^ K in ~ 10® yr. This is because the pressure 
rise makes the destruction time of the majority of molecular clouds in the galaxy 
less than 10® yr. The tidal force from the potential well of the cluster accelerates 
molecular clouds in a disk galaxy infalling towards the cluster center. Thus, 
the kinetic pressure rises above ~ 3 x 10^ cm~^ K. Before the galaxy reaches 
the cluster center, the star formation rate reaches a maximum. The peak is 
three to four times larger than the initial value. If this is the main mechanism 
of the Butcher-Oemler effect, blue galaxies are expected to be located within 
~ 300 kpc from the center of the cluster. However this prediction is inconsistent 
with the recent observations. The increase of external pressure when the galaxy 
plows into the intracluster medium does not change star formation rate of a 
disk galaxy significantly. Thus, the increase of external pressure may not be 
the main mechanism of the Butcher-Oemler effect. The velocity perturbation 
induced by a single high-speed encounter between galaxies is too small to affect 
star formation rate of a disk galaxy, while successive high-speed encounters 
(galaxy harassment) trigger star formation activity because of the accumulation 
of gas in the galaxy center. Therefore, the galaxy harassment remains as the 
candidate for a mechanism of the Butcher-Oemler effect. 
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1. Introduction 

Clusters of galaxies in the redshift range of 0.2 — 0.5 often exhibit an overabundance, 
relative to present-day clusters, of blue galaxies (Butcher & Oemler 1978). This star 
formation activity is often called the Butcher-Oemler effect (BOE). Several mechanisms 
responsible for the effect have been suggested. 

Tidal compression by the cluster potential increases the velocity dispersion of clouds 
in a disk galaxy (Byrd & Valtonen 1990 ; Valluri 1993). Henriksen & Byrd (1996) estimate 
the influence of tidal acceleration by the cluster potential on galaxies. They show that if 
the cluster potential is centrally peaked, compression of galaxy disk is significant. They 
find that the velocity perturbation due to the gravitational tidal forces are large enough for 
molecular clouds to give birth to stars. 

Dressier & Gunn (1983) suggest that a significant fraction of galaxies in a cluster 
undergo a rapid rise in external pressure as they plow into the hot intracluster medium. 
Evrard (1991) numerically confirms this possibihty. He suggests that the high-pressure 
leads to star formation activity. 

Recently, Moore et al. (1996) numerically show that successive high speed encounters, 
"galaxy harassment", cause a complete morphological transformation from disks to 
spheroidals. Since galaxy harassment is an extremely effective mechanism for fueling central 
region of a galaxy, it may cause star formation activity. 

However, the above studies do not derive the history of star formation rate (SFR) 
quantitatively. In this paper, a simple method is devised to evaluate the SFR of a disk 
galaxy in a cluster in order to find out what mechanism drives the BOE. This model is 
useful to predict color and spatial distributions of the blue population in BOE clusters. 
The above three mechanisms, (1) tidal force from the potential well of the cluster, (2) 
increase of external pressure, and (3) high-speed encounters between galaxies are studied as 
candidates of the BOE. These mechanisms provide pressure variations of local interstellar 
medium (ISM) in a disk galaxy because of increase of kinetic or thermal pressure of ISM, 
and gas inflow towards the galaxy center. As will be shown in detail, the time scales of 
the variations are ~ 10^, which reflect the dynamical time scales of a galaxy and central 
region of a cluster. Star formation responds this pressure variation locally. Elmegreen & 
Efremov (1997) obtain a simple relation between star formation efficiency of molecular 
clouds of varying mass and pressure. Their results indicate that if the pressure of a cloud 
increases the order of a magnitude, the star formation efficiency increases and the lifetime 
decreases significantly. The relation does not care about the origin of the pressure. Thus, 
it can be applied to the above three cases (1) - (3). Since the lifetime of clouds is ^ 10® yr 
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(Elmegreen & Efremov 1997), the pressure variations in the cases (1) - (3) are expected to 
affect the evolution of all the molecular clouds and the SFR of a disk galaxy. This will be 
discussed in the following sections. 

The plan of this paper is as follows. In §11, a model of molecular cloud evolution is 
described. In §111, the pressure evolution of a disk galaxy in a cluster is derived, and the 
SFR is evaluated. Finally, §IV is devoted to conclusions. In that section, observational 
implications are commented. 



2. Molecular Cloud Models 

Molecular clouds are divided according to their initial masses, that is, 
-^min = Ml < ... < Mi < ... < Mmax, whcrc the lower and upper cutoffs are 

-^min = 10^ Mq and Mmax = 10^'^ ^0) respectively, and their intervals are chosen that 
log(M(+i/Mj) = 0.01. Referring to the total mass of clouds whose initial masses are between 
Mi and Mj+i as AMj(t), the rate of change is 

'-^--MS.,HS^M-^, (1) 

where /j is the initial mass fraction of the molecular clouds whose initial masses are between 
Mi and Mj+i, is the gas ejection rate from stars, 5'inoi(^) is the recycle rate of molecular 
gas, T{Mi,P) is the destruction time of a molecular cloud with mass Mi and pressure P, 
respectively. 

Since the mass spectra of open clusters, globular clusters, and molecular clouds follow 
N{M) (X M-"", where a = 1.6 - 2 (Harris & Pudritz 1994 ; Elmegreen & Efremov 1997), 
the initial mass function of molecular clouds is assumed to be N{M) oc M"^. Thus fi can 
be derived from this relation. 

The gas ejection rate from stars is divided into two terms, S'*(t) = Ss{t) + Si. The gas 
ejection rate from stars with small lifetime is given by 

rrriu 

Ss{t)^ ip{t-tm)R{m)(l){m)dm, (2) 

J mi 

where m is the stellar mass, V^(t) is the SFR, 0(m) is the IMF expressed in the form of the 
mass fraction, R{m) is the return mass fraction, and tm is the lifetime of stars with mass 
m. The slope of the IMF is taken to be 1.35 (Salpeter mass function). The upper and 
lower mass limits, and m/, are taken to be 50Mq and 2.8Mq, respectively, although the 
lower limit of the IMF is O.OSMq. The hfetime of stars and the return mass fraction are 
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adopted from those in Mihara & Takahara (1994) and Maeder (1992), respectively. Note 
that the hfetime of a star with 2.8Mq is 4 x 10® yr, which is identical with the end time of 
calculations in the later sections. The gas ejection rate from stars with large lifetime, Si, 
will be specified later. 

The recycle rate of molecular gas is 

Sr.o^{t) = U^-e{M,,P)]^^^ (3) 

where e(Mj, P) is the star formation efficiency of a molecular cloud with mass Mj and 
pressure P. Under simple assumptions, Elmegreen & Efremov (1997) derive e{Mi,P) and 
T{Mi,P) : their results are adopted in the following. Note that the pressure P includes the 
kinetic pressure (Elmegreen 1989 ; Elmegreen & Efremov 1997). In our Galaxy, the kinetic 
pressure is larger than thermal pressure (Elmegreen 1989 ; Elmegreen & Efremov 1997). 



The SFR is described by 



m = E<M'.n^^. (4) 



Equation ([1|) can be solved if P is given. In the next section, P is determined for 
several cases. Then the SFRs are calculated form Equations (0) — (|). 



3. Evolution of SFR of Disk Galaxies in a Clusters of Galaxies 

3.1. General Analysis 

In this subsection, the evolution of pressure is treated as a parameter to find the 
relation between the evolution of pressure and the SFR of a disk galaxy. 

Three cases are considered, that is, the final value of P is Pf = 10, 100, and 1000 Pq 
(Figure 1(a)). The pressure is given in units of Pq = 3 x 10^ cm~^ K which is the total 
pressure in the solar neighborhood (Elmegreen & Efremov 1997). The pressure P increases 
for t > 0. For t < 0, the model galaxy is assumed to be at P = IP© 

The calculation starts at t = —10^ yr. At that time, the total mass of molecular clouds 
in the model galaxy, Mtot, is 10^'^ Mq, and the mass spectrum of molecular clouds is 
cx M~^. For t < 0, ignoring Equation (^, the source term S"* in Equation (Q) is artificially 
fixed at 5*0 = 6 Mq yr^^ which is the SFR of our Galaxy (Gusten & Mezger 1983). After 
the calculation starts, the mass spectrum is modified because the destruction time of 
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molecular clouds depends on their mass (Elmegreen & Efremov 1997). Until t — 0, the 

destruction and formation rates of molecular clouds become identical to within 3%. At 
t = 0, Mtot = 2.4 X 10^ ^^^d the mass spectrum of molecular clouds is cx M~^-^. The 
total mass of molecular clouds is similar to that of our Galaxy (~ 2 x 10^ Mq, Larson 
1987). Note that for t > the SFR history normalized by 5*0 is independent of the value of 
5*0. The gas ejection rate from stars with large lifetime is given by Si = So — Ss{0) for t > 
and for fixed Pq. 

The evolution of pressure, SFR, and the total mass of molecular clouds arc indicated 
in Figure 1. As is shown, the SFRs and the time at which it become maximum depends 
on P{t). Taking the case where Pf = 10 Pq as an example, the SFR can be explained 
as follows. The pressure P reaches 10 Pq at t — 10^ yr. When P = 10 Pq, the mass 
of a molecular cloud whose destruction time is less than 10^ yr is ^ 10^ Mq (see Figure 
4 in Elmegreen & Efremov 1997). On the other hand, when t — 0, the total mass of 
molecular clouds whose individual mass is less than 10^ Mq is 1.6 x 10^ Mq. The ratio, 
1.6 X 10^ Mq/10^ yr = 16 Mq yr~^, nearly corresponds to the peak of the SFR (Figure 
1(b)). In the case where Pf = 100 Pq, the SFR reaches a maximum before P becomes 
100 Pq (t = 2 X 10* yr). This is because all molecular clouds (< 10*^-^ Mq) are affected 
by the pressure increase before t = 2 x 10* yr, since the destruction time of a molecular 
cloud with M = 10*-^ Mq is ~ 10* yr when P = 100 Pq. In fact, Mtot is significantly 
reduced when t — 2 x 10* (Figure 1(c)). Therefore, the peak of the SFR is determined by 
the balance between the time scale of pressure increase and the destruction time of clouds. 
Moreover, Figure 1(b) shows that the evolution of SFR is not very sensitive to P/, if it is 
larger than 100 Pq and if the time scale of pressure increase is ~ 10* yr. This is because 
almost all of the clouds are destroyed regardless of Pf. 

In summary, when pressure increases on a time scale of 10* yr, the SFR rises 
significantly if P/ ^ 10 Pq but reaches the ceiling at P/ ~ 100 Pq. The above arguments 
should be true for all of the BOE mechanisms treated in the following sections, because this 
model does not care about the origin of pressure. 

3.2. A Radially Infalling Galaxy 

A radially infaUing disk galaxy is affected by tidal force from the cluster potential well. 
Molecular clouds in the galaxy are accelerated. Thus the pressure P from random cloud 
motions increases. 
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Three types of mass profiles in the central region of the cluster are investigated; 

/ ^2 \ 1/21 



MrA< r) = 1.2 X 10^^ M, 
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Equations (^) is the profile inferred from gravitational lensing constraints, Equations 
is the so-called King model profile, and Equations (0) is the profile obtained in numerical 
simulations done by Navarro, Frenk, & White (1995, 1996). The characteristic radii for 
each are tgl = 50 kpc, rxing = 200 kpc, and tnew = 900 kpc, respectively. The first profile 
is the same as that in Henriksen & Byrd (1996). The mass distributions are normalized so 
that the mass in side 500 kpc is 2.5 x 10^^ Mq for all distributions. When the mass profile 
of the cluster follows Equation (|^, the profile of the intracluster gas can be fitted well by 
(^) except for the normalization of mass. In that case, the ratio of characteristic radii, 
'^Kmg/''"NFW5 is 0.22 (Makino, Sasaki, & Suto 1998), which is the reason for taking a radius 
of 900 kpc in Equation (^. 

In order to calculate the tidal acceleration and the velocity perturbation of molecular 
clouds, the impulse approximation is used as in Henriksen & Byrd (1996). Strictly speaking, 
this approximation may not be good because the rotation period of a galaxy is comparable 
to the time-scale in which the galaxy crosses the central region of the cluster. However, as 
will be noted, the results are not very sensitive to the velocity perturbations for the models 
considered in this subsection. 

The transverse and radial tidal accelerations are 

R 

at = GMJ< r) 



GMJ<r) 



(r + i?)2]3/2 

1 

' (r + Rf 



(9) 



where G is the gravitational constant and R is the radius of the galaxy (= 20 kpc). The 
velocity perturbation can be estimated as 

V= t adt , (10) 



where a is the transverse or radial tidal acceleration. The kinetic pressure is given by 
-Ptidai = PdV"^, where the average density of the galaxy disk is pd = 1-67 x 10^24 g cni-3. 
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The evolution of a is calculated from the orbits of the galaxy in the potential well of the 
cluster. As an idealized case, a radially infalling galaxy is investigated. The velocity of the 
galaxy, t>gai, can be obtained by solving the equation of motion. At t = 0, the position of 
the galaxy is 350 kpc away from the center of the cluster. From there the galaxy radially 
infalls towards the center. The initial velocity is 1000 km s~^. Figure 2 shows the orbits of 
the galaxy. 

Figures 3(a) and 4(a) show the evolution of the pressure, P = Pq + Ptidai, for the three 
types of potentials (Equations |^, and The rates of pressure increase due to the radial 
tidal acceleration are somewhat larger than those due to the transverse tidal acceleration. 
Since the galaxy rotates, and radial and transverse tidal acceleration cancel, the pressure 
of a real galaxy is less than that shown in Figure 4(a). Note that when the galaxy disk is 
perpendicular to the orbit. Figure 3(a) indicates the real tidal acceleration because there is 
no radial tidal acceleration. 

Using the pressure evolution. Equation (|I]) is solved. The evolution of a model galaxy 
for t < is the same as that in § |3.1| . The evolution of the SFRs is shown in Figure 3(b) 
and 4(b). The SFRs become maximum before the galaxy reaches the center of the cluster. 
At that time, the galaxy is located at 50 ^ r 200 kpc (Figure 2). Since the time scales of 
the pressure increase are ~ 2 x 10*^ yr and P becomes larger then 100 Pq for all models. 



the evolution of the SFRs is similar among them (see § |3.1| ). Thus, it is common among the 



models that the maximum SFR is three to four times as much as the initial SFR. 

The evolution of Mtot is shown in Figure 3(c) and 4(c). During the passage through 
the central region of a cluster, Mtot decreases to ^ 1/5 of the initial value. 



3.3. Pressure from the Intracluster Medium 

Numerical simulations of rich cluster evolution show that the local static pressure 
around a significant fraction of galaxies experienced a rapid increase at high redshift (Evrard 
1991). Thus, the change of SFR due to this pressure increase is investigated. The evolution 
of the pressure is shown in Figure 5(a) ; for < t < 3 X 10* yr, P = 10^*/^^ Pq, and for 
t > 3 X 10* yr, P = 10°'^ Pq. This corresponds to the most rapid increase of pressure in 
Figure 1 of Evrard (1991). The evolution of a model galaxy for t < is the same as that 



in § |3.1| . Figure 5 also shows the evolution of SFR and Mtot- The SFR increases by at most 



70%. Thus, the increase of static pressure may not induce observable changes in color. 
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3.4. High-Speed Encounters 

Moore et al. (1996) propose that muhiple high speed encounters between galaxies 
(galaxy harassment) change the structure of galaxies in clusters. Therefore, the effect of 
high-speed encounters on the SFR of a disk galaxy is considered. 

In the first place, a single encounter is considered. As an extreme penetrating 



encounter is calculated. The derivation of P{t) is the same as that in § |3.2| except for the 
mass distribution of a perturber. 



The mass distribution of a perturber is 

„2 
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for r < 50 kpc. The characteristic radius rper is 5 kpc. For r > 50 kpc, 
Mper(< t) = 2.0 X 10^^ Mq. The evolution of a model galaxy for t < is also the 
same as that in § |3.1| . At t = 0, a model disk galaxy is located at r = 100 kpc, and the 
relative velocity to the perturber is 1000 km s~^. Because of the symmetry between the 
approach and recession phases of the model galaxy, only the transverse tidal acceleration is 
calculated. The results are shown in Figure 5. The model galaxy passes the center of the 
perturber at t ~ 10^ yr. Since the kinetic pressure does not change significantly, the SFR 
rises by only 13%. 

In this argument, however, the ISM of the perturber is not taken into account. If it has 
its own ISM with density pp, the ram pressure ppV"^ acts on molecular clouds in the model 
galaxy. For example, when = 1.67 x 10"^'' g cm~'^ for r < 5 kpc and pd = for r > 5 
kpc, the SFR temporarily increases more than ten times larger than the initial value. The 
formation of a ring or tidal arms due to a close encounter will also affect the SFR. However, 
Moore et al. (1996) indicate that such close encounters (< 30 kpc) are extremely rare in a 
cluster. Thus, a single high speed encounter does not drive the BOE. 

On the other hand, galaxy harassment would induce star formation activity in many 
galaxies in a cluster. Using the results of Moore, Lake, and Katz (1998) and Lake, Katz, & 
Moore (1998), the change of SFR due to the harassment is calculated. 

Galaxy harassment has a great influence on galaxies smaller than the Galaxy. Thus, 
the initial conditions of a model galaxy are different from those considered above. At t = 0, 
the SFR is 1.2 Mq yr~^ and Mtot = 3.7 x 10^ Mq, which correspond to those of the model 
galaxy in Moore et al. (1998) and Lake et al. (1998). The initial value of pressure is taken 
to be P = 0.1 Pq in order to balance the cloud destruction rate with the formation rate. 
Lake et al. (1998) show that most of the gas of the harassed galaxy is driven to the region 
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within 1 kpc of the center within 3 Gyr. In the extreme case, half of the mass is transferee! 
in an interval of 100 - 200 Myr. Assuming that 60% of the gas is driven to the region 
within 1 kpc of the center and that the hight of the gas disk is 200 kpc, the gas density is 
2.4 X 10~^^ g cm~^. Since the effective velocity is ~ 130 km s~^ (Fig. 2 in Moore at al. 
[1998]), the pressure of the central region of the harassed galaxy is 1.0 x 10"^ Pq. Thus, the 
evolution of the model galaxy is assumed to be P{t) = IQ^^^^/'^y^)'^ Pq (Figure 5a). Using 
this and Equation (|I|), the evolution of the model galaxy can be calculated. The results 
are shown in Figure 5(b) and 5(c). The history of SFR shows that galaxy harassment can 
induce strong star formation activity. 

4. Conclusions 

In this paper, the influence of environmental effects on the star formation rate (SFR) 
of disk galaxies in clusters was quantitatively estimated. The constructed model treats the 
response of star formation to local pressure variance. Since it does not care about the origin 
of pressure, it can be used to investigate various candidates for the Butcher-Oemler effect. 

General analysis indicates that the SFR increases significantly when the pressure of 
molecular clouds rises above ~ 10 Pq in ~ 10^ yr but reaches the ceiling at ~ 100 Pq. 
This is because the destruction time of most molecular clouds in the galaxy becomes ^ 10® 
yr. 

In the analysis of a disk galaxy tidally affected by the cluster, it is found that the 
maximum SFR is around three or four times as much as that of a field galaxy, and does 
not depend greatly on the mass profiles of the cluster and the type of the tidal force. The 
SFR becomes maximum before the galaxy reaches the center of the cluster. At that time, 
the distance from the center is 50 ^ r 200. Fritze -v Alvensleben & Gerhard (1994) have 
investigated the color evolution during and after starbursts for many types of galaxy. Their 
results show that the color of a galaxy significantly becomes bluer for ~ 10® yr after the 
burst if the gas mass fraction just before the burst is ~ 10 %, the duration of the burst 
is ~ 10® yr, and most of the gas is consumed during the burst. Therefore, if the tidally 
affected galaxies are responsible for the Butcher-Oemler effect (BOE), the blue galaxies 
should be found ^ 300 kpc from the center of the cluster (Figure 2). It is confirmed that 
even if the initial velocity of the model galaxy is 2000 km s^^, which is larger than the virial 
velocity of the model cluster, this prediction does not significantly change. However, recent 
observations are inconsistent with this prediction. In A2317, the blue population decreases 
remarkably inside 300 kpc (Rakos, Odell, & Schombert 1997). Since this trend is observed 
in other distant clusters (Abraham et al. 1996 ; Balogh et al. 1997), it appears that tidal 
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force from the potential well of a cluster does not induce the BOE. 

As a galaxy plows into the hot intracluster medium, the external pressure rapidly rises. 
However, the rise of local pressure based on numerical calculation done by Evrard (1991) 
does not appear to induce starbursts because the maximum of the pressure is small. This is 
consistent with the recent observation of a merging cluster (Tomita et al. 1996). 

The velocity perturbation induced by a single high-speed encounter is too small 
to affect the SFR of a disk galaxy. However, successive high-speed encounters between 
galaxies (galaxy harassment) lead to gas inflow and strong star formation activity. Thus, 
galaxy harassment may be a likely explanation for the star-forming galaxies observed in 
high-rcdshift clusters (Oemler, Dressier, & Butcher 1997 ; Couch et al. 1998 ; Small et 
al. 1998). In order to know whether galaxy harassment is the m,am origin of BOE, it is 
required to estimate the fraction of galaxies in which the gas accumulates in the center. 
The spatial distribution of blue galaxies predicted by the galaxy harassment model may be 
consistent with the observations, because galaxies in the core of clusters will be older than 
galaxies at the edges. 

On the other hand, the BOE may not be excess star formation induced by the infall 
process or internal tides in clusters of galaxies. Abraham et al. (1996) suggest that the 
BOE is due to the increased rate of infall of bluer field galaxies at higher redshift, and 
that the star formation is truncated without an increase. A model including ram-pressure 
stripping and evolution of intracluster medium will be needed to investigate this idea. 

I would hke to thank S. Sasaki and T. Tsuchiya for their helpful comments. I wish to 
thank S. Inoue for his assistance. I am also grateful to an anonymous referee for several 
suggestions that improved this paper. This work was supported in part by the JSPS 
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Figure Captions 

Fig.l — Histories of (a) pressure (b) star formation rate, and (c) total mass of 
molecular clouds. 

Fig. 2 — - Orbits of radially infalling galaxies in the three potential wells. Solid line; 
lensing potential (GL), dotted line; King model (King), dashed line; the one derived by 
Navarro et al.(NFW). 

Fig. 3 — - Histories of (a) pressure (b) star formation rate, and (c) total mass of 
molecular clouds for a galaxy affected by the transverse tidal acceleration. Solid line; lensing 
potential (GL), dotted line; King model (King), dashed line; the one derived by Navarro et 
al.(NFW). 

Fig.4 — Same as in Fig.3 but for radial tidal acceleration. 

Fig.5 — Histories of (a) pressure (b) star formation rate, and (c) total mass of 
molecular clouds for a galaxy affected by the pressure of intracluster medium (ICM; sohd 
line), a single high-speed encounter (HE; dotted hue), and galaxy harassment (GH; dashed 
line) . 
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